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We present a systematic study of the ballistic electron conductance through sp and 3d transition 
metal atoms attached to copper and palladium crystalline electrodes. We employ the ab initio 
screened Korringa-Kohn-Rostoker Green's function method to calculate the electronic structure 
of nanocontacts while the ballistic transmission and conductance eigenchannels were obtained by 
means of the Kubo approach as formulated by Baranger and Stone. We demonstrate that the 
conductance of the systems is mainly determined by the electronic properties of the atom bridging 
the macroscopic leads. We classify the conducting eigenchannels according to the atomic orbitals of 
the contact atom and the irreducible representations of the symmetry point group of the system that 
leads to the microscopic understanding of the conductance. We show that if impurity resonances 
in the density of states of the contact atom appear at the Fermi energy, additional channels of 
appropriate symmetry could open. On the other hand the transmission of the existing channels 
could be blocked by impurity scattering. 

PACS numbers; 73.63.Rt, 73.23. Ad, 75.47.Jn, 73.40.Cg 



I. INTRODUCTION 

State-of-the-art experimental techniques, like scanning 
tunneling microscopy (STM)- and mechanically control- 
lable break junctions (MCBJ),^ "^ enable the manipula- 
tion of individual atoms and the fabrication of metal- 
lic, atomic-sized contacts. From a fundamental scientific 
point of view these systems allow us to investigate a va- 
riety of phenomena at the ultimate limit of atomic scale 
and prove basic ideas of quantum mechanics. For ex- 
ample, when studying conductance, the familiar Ohm's 
law completely breaks down and the fully quantum- 
mechanical description should be applied to elucidate the 
phenomenon j'^'^1^1^1^1^ 
I Our understanding of ballistic transport of atomic- 
sized conductors is usually based on the Landauer 
formulafiSi 




Tl=l 



where conductance is represented as a sum of transmis- 
sion probabilities r„ of individual eigenchannels. When 
the TnS are known many properties of a system can 
be predicted such as conductance fluctuationsfii the 
shot noiseji^ dynamical Coulomb blockadai^ and the 
supercurrent^ It has been demonstrated in the pio- 
neering work by Scheer et al^ that a study of the 
current-voltage characteristics for the superconducting 
atomic-sized contacts allows to obtain transmission prob- 
abilities T„'s for particular atomic configurations real- 
ized in the MCBJ experiments. The r„'s can be de- 
termined by fitting theoretical and experimental I — V 
curves and exploiting information contained in the " sub- 
gap structure" >iM6 The breakthrough in the interpreta- 
tion of experimental results was done by Cuevas et alA 
who suggested a parameterized tight-binding model for 



the description of atomic constrictions. As a central re- 
sult it has been shown that the number of conducting 
modes (eigenchannels) for a single-atomic contact cor- 
relates with the number of valence orbitals of the con- 
tact atom which therefore determines the current fiow- 
ing through the system. This result has been confirmed 
experimentally by Scheer et alm^ The analysis of the last 
conductance plateau corresponding to single-atom con- 
strictions has revealed that a number of conducting chan- 
nels is 1 for Au, 3 for Al and Pb, and 5 for Nb. 

The main goal of this paper is to clarify with the help of 
parameter-free ab initio calculations the issue on how the 
valence states of the central atom of a single-atom contact 
determine conducting channels. Our approach is based 
on the screened Korringa-Kohn-Rostoker (KKR) Green's 
function methodii supplemented with the Baranger and 
Stonc^* formulation of the conductance problem. In- 
stead of considering a variety of different metallic atomic 
point contacts we restrict our calculations to single sp 
or 3d transition metal atoms attached to noble metal 
(Cu) or transition metal (Pd) leads. Systems of this 
type can be realized experimentally by using alloys to 
fabricate the point contacts. ^^'^"^ According to our pre- 
liminary studyjSi we await the variation of conductance 
depending on the atomic number of the impurity situated 
at the central site. The number of conducting modes is 
also expected to change with respect to the atomic num- 
ber of the impurity, especially because of the existence of 
the virtual bound states. We would like to mention, that 
the conductance through different type of constrictions 
contaminated by impurities has been studied before on 
the ab initio level by Lang^, Hirose and coworkers^'' and 
Palotas and coworkersi24 In contrast to previous stud- 
ies, we analyze the eigenchannel decomposition of the 
conductance. Within our approach each channel is clas- 
sified according to the irreducible representation of the 
corresponding point group of the system as well as by 



orbital contributions when the channel wave function is 
projected on the central atom. That relates the method 
presented here with the tight-binding model suggested 
by Cuevas et al& and gives a microscopic insight into the 
conductance. 



II. THEORETICAL APPROACH 

A. Set-up of the problem 

The structures under study consist of two semi-infinite 
Cu or Pd fee (001) leads connected by an atomic cluster 
with a single contact atom in the middle of a junction as 
it is shown in Fig.l. We consider the contact atom to be 
either a host atom or different sp or 3d transition metal 
impurities as will be discussed later. We assume that the 
geometry of constrictions shown in Fig. 1 catches limiting 
configurations of atomic contacts produced in the MCBJ 
experiments. To realize different contact atoms one could 
use dilute alloys as a wire material in the MCBJ^iSiSS or 
probe conductance through a single atom with an STM 
tip.^^ Since in the experiments the measurements are 
not performed for the equilibrium atomic configurations 
rather than the conductance is recorded dynamically un- 
der stretching of nanocontactSyS- we did not optimize the 
interatomic distances and kept them the same as in bulk 
fee metals Cu and Pd with lattice constants 6.83 a.u. and 
7.35 a.u., respectively. 



B. Electronic structure calculation of atomic 
contacts 

Our calculations are based on density functional theory 
within the local density approximation. We employed the 
non-relativistic version of the screened Korringa-Kohn- 
Rostoker (KKR) Green's function method to calculate 
the electronic structure of the systems. Since details of 
the method can be found elsewheret^S, only a brief de- 
scription is given below. The parametrization of Vosko, 
Wilk, and Nusair^^ for the exchange and correlation en- 
ergy was used. The potentials were assumed to be spher- 
ically symmetric around each atom (atomic sphere ap- 
proximation, ASA). However, the full charge density, 
rather than its spherically symmetric part, was taken into 
account. To achieve well converged results the angular 
momentum cut-off for the wavefunctions and the Green's 
function was chosen to be ^max — 3 that imposed a nat- 
ural cut-off 2Zniax = 6 for the charge density expansion. 
In case of the heavy element Pd the scalar relativistic 
approximation-- was used. 

In the multiple-scattering KKR approach the one- 
electron retarded Green's function is given by the site 




FIG. 1: Geometry of Cu and Pd atomic-sized constrictions 
studied in the paper: (a) three atom chain suspended between 
the leads; (b) pyramid- like contact. The central (contact) 
atom of the constriction can be either a host or impurity atom. 
The conductance is calculated between left (Sl) and right 
(Sk) planes positioned in the leads. 



angular momentum expansion 

G+(R„+ r,R„, + r';E) 
= 6nn'^Y.^L{r<;E)H2{r>;E) (1) 

L 

+ Y.Rl{v-E)Gli,{E)Ri,{v'-E) 

LL' 

where r, r' are restricted to the cells n and n'; r<, r> 
denote one of the two vectors r or r' with the smaller or 
the larger absolute value, and local functions R^{y\E) 
and H2{v]E) are the regular and irregular solutions of 
the Schrodinger equation for the single potential y,i(r) 
of the rt-th cell in free space. Here the index L — {I, m) 
stands for the angular momentum quantum numbers and 
atomic units are used: e = — \/2, h — 1, m = 1/2. 
The structural Green's function G'll, (E) (structure con- 
stants) in Eq. ^ is related to the known structure con- 
stants of the appropriately chosen reference system by 
the algebraic Dyson equation which includes the differ- 
ence At = Snn'^LL' between local ^-matrices of the 
physical and a reference system. In the screened KKR 
methodic we use a lattice of strongly repulsive, constant 
muffin-tin potentials (typically, ~ 4Ry height) as refer- 
ence system that leads to structure constants which decay 
exponentially in real space. 

When the KKR method is applied to the systems 
shown in Fig.l both the constriction region and the leads 
are treated on the same footing. This is achieved by us- 
ing the hierarchy of Green's functions connected by a 
Dyson equation, so that one performs the self-consistent 
electronic structure calculations of complicated systems 
in a step-like manner. First, using the concept of princi- 
pal layers together with the decimation technique)^ we 
calculate the structural Green's function G^jf/ (k|| , E) of 
the auxiliary system consisting of semi-infinite leads sep- 
arated by a vacuum barrier (here j and j ' refer to layer 
indexes) . At the the second step the atomic cluster is em- 
bedded between the leads by solving the Dyson equation 



3 



self-consistently, 



D. Conductance eigenchannels 
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and the structural Green's function GJ^' {E) of the com- 
plete system is obtained. Eq. (0) is solved in real space 
since due to effective screening of the perturbation the 
charge deviations are restricted nearby the constriction. 



C. Evaluation of conductance 

In what follows, wc will consider only a ballistic co- 
herent regime of the electron transport thus limiting 
our study to zero temperature, infinitely small bias and 
singe-particle picture. To evaluate the conductance of 
nanocontacts we employ the Kubo linear response the- 
ory as formulated by Baranger and Stone;^^. 

g = Go f dS I dS' 

X G+ir,r',EF) d,d.' G-{v\v,Ef), (3) 

here Gq — 2e^//i, G^ and G+ are retarded and advanced 
Green's functions, respectively. The current flows in z di- 
rection, and f dz g = f{dzg) — {dzf)g- The integration 
is performed over two (left and right) planes S'l and S'r 
which connect the leads with the scattering region as it 
is shown in Fig. 1 . The implementation of the Eq. JJl in 
the site angular momentum representation of the KKR 
method and related to it convergence properties were de- 
scribed in detail in Ref. Within the ASA used in 
this work instead of integration over planes we average 
conductance over atomic layers and obtain: 



g 



GoTr(„,i) [Di^GD^G^] 



(4) 



Here, Tr involves site (n) and angular momentum [L) 
indices related to the atomic plane S'l (Fig. 1). G = 
{GYl {^f)} stands for the matrix notation of the struc- 
tural Green's function introduced in Eq. (Q, taken at 
the Fermi energy. Site-diagonal operators Di^ and I?r 
related to the left (L) and right (R) leads are defined as 



l(R)Jll' 
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,(5) 



where different sings, and "— ", refer to different op- 
erators, Dl and Dr, respectively. The integral is per- 
formed over the volume V^^^^-^ of the Wigner-Seitz cell 
around site n which belongs either to the atomic plane 
S'l or to the plane Sr, and A is the distance separating 
nearest-neighbor atomic layers in z direction. 



Ballistic conductance of an atomic constriction can 
be decomposed into individual eigenchannels. Accord- 
ing to the Landauer formalism, which was proved to be 
equivalent to the Kubo approachfiS conductance reads 
a.s g — GoJ2n'''ni-^F)- The r„'s are eigenvalues of the 
matrix T'k/i:k>' which defines the transmission probabil- 
ity for an incident Bloch wave from the left electrode to 
be transmitted to an outgoing Bloch wave in the right, 
where k and k' are Bloch vectors for incident and trans- 
mitted waves, and /i and /i' are the corresponding band 
indices. A procedure for the evaluation of the conduc- 
tance eigenchannels will be published elsewhere. '^^ To 
solve the eigenvalue problem for the transmission ma- 
trix T, the perturbed Bloch states of the whole system 
are projected to the basis of local functions BJl{v) and 
conductance is represented as a sum of two contributions: 



g = GoTr(,,i) [T+] + Go Tr [T"] 



with 



r± - ± J±Dt GDnGU ±Dt 



(6) 



(7) 



Here are two parts of the anti-symmetric hermitian 
operator Dl = + defined by Eq.©. The and 
Dj^ have positive and negative eigenvalues, respectively. 

When the atomic plane S'l is placed in the asymptotic 
region of the left lead far from the atomic constriction the 
second term in Eq. © , involving trace of T~ , is equal to 
zero and eigenchannels are found as eigenvectors of op- 
erator T+. However, in practice we are forced to take 
integration planes closer to the constriction in order to 
obtain convergent value for the conductance with respect 
to number of atoms in the cross-sections of the leads. The 
positions of the atomic planes Sl and Sr, therefore, do 
not meet the asymptotic limit criterium as it is formally 
required by the scattering approach. However, since the 
current through the structure is conserved, any position 
of the planes is suitable for the calculation of conduc- 
tance. Thus, if Sl is placed somewhere in the scattering 
region we have to sum up all multiple scattering contribu- 
tions. All contributions in direction of the current cause 
T"^, whereas all multiple scattering contributions in op- 
posite direction of the current give rise to T~. In the 
region of the lead where the potential is a small pertur- 
bation with respect to the bulk potential the contribution 
to the conductance due to T~ is one order of magnitude 
smaller than T+. To identify the transmission proba- 
bilities of individual eigenchannels the spectrums of the 
operators T"*" and T~ are arranged in a proper way us- 
ing the symmetry analysis of eigenvectors. For further 
details we refer to Ref. l3ll 

Most of the structures studied in this work obey 
the symmetry determined by the fee (001) surface 
(Fig.l). Further we denote individual channels by the 
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FIG. 3: Symmetry resolved LDOS at impurity atoms placed 
at the contact site of Cu junction [see Fig. 1(a)]. 



FIG. 2: Top panel (a): conductance and corresponding eigen- 
channel decomposition through the atomic chains with sp im- 
purities connected to the Cu leads [see Fig. 1(a)]. Bottom 
panel (b): Symmetry resolved LDOS at impurity sites, at the 
Fermi energy, vs atomic number of impurity. 



indices of the irreducible representations of this group us- 
ing notations of Ref. 133 . common in the band theory. In 
addition, each channel can be classified according to the 
angular momentum contributions when the wave func- 
tion of the channel is projected on the contact atom of 
the constriction. This is very helpful since the channel 
transmission can be related to the states of the contact 
atomi^ For example, the identity representation Ai of the 
Civ group is compatible with the s, Pz and d^^ orbitals 
(here z is the axis perpendicular to the surface), while 
the the two-dimensional representation A5 is compatible 
with the Px, Py, dxz, dyz orbitals. 



III. RESULTS AND DISCUSSION 

A. Conductance through sp atoms 

We consider first atomic constrictions shown in Fig. la 
which are simulated by the straight three atom chains 
suspended between copper leads with different impuri- 
ties in the center of the junction. The calculated conduc- 
tance value of a pure Cu nanocontact is 0.96 Gq. Only 
one channel contributes to the total transmission, it has 
Al symmetry and is projected on s and Pz orbitals of the 



contact Cu atom. In addition, small contributions from 
states with higher angular momentum {dz2 and /^a) are 
also present. This result is in agreement with other the- 
oretical studies^ as well as with experiments for noble 
metalsiSiM 

We proceed with the discussion of the conducting prop- 
erties of single sp atoms {Z ^ 11 . . . 16) attached to Cu 
electrodes (Fig. la). In Fig. 2 (top panel) we present the 
conductance together with the eigenchannel decomposi- 
tion and compare it with the symmetry projected local 
density of states (LDOS) at the impurity atom, at the 
Fermi energy (bottom panel). We observe from Fig. 2a 
that the conductance changes significantly depending on 
the type of sp impurity bridging the Cu leads. For ex- 
ample, conductance through the sodium atom is smaller 
than 0.5 Gq while in case of phosphorus it reaches a value 
of 2.4 00. Our calculations show that three channels are 
involved in conductance: one sp^-channel of Ai symme- 
try and a double degenerate channel of A5 symmetry 
which locally consists of Px, Py states when projection 
on the contact atom is performed. From Fig. 2 we see 
that the variation of the eigenchannels transmission with 
atomic number follows the changes in the orbital resolved 
LDOS at Ef- 

In case of Na, Mg and Al the transmission is essen- 
tially determined by one highly symmetric Ai channel 
(Fig. 2a). For sodium the electronic states at Ep have 
mainly s character (Fig. 3) but transmission is small be- 
cause of weak coupling to the Cu leads. As we move up to 
Mg and Al the s states are filled up with more electrons. 
At the same time, the Pz states contribute significantly to 
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FIG. 4: Conductance and its eigenchannel decomposition for 
the ziz-zag atomic chains suspended between the Cu leads 
with sp atoms placed at the contact site. 



the LDOS at Ep (Fig. 3). That causes a saturation of the 

Al channel transmission. For larger atomic numbers the 
Al channel is determined by states only. The trans- 
mission of the Al channc^l is reduced again since for Si, 
P and S impurities the s states are occupied. The results 
of the calculations shown in Figs. 2 and 3 reveal that the 
increase of conductance through Si, P and S impurities is 
due to a virtual bound state formed by Px and Py orbitals 
(dashed line in Fig. 3). This state gives rise to a double 
degenerate A5 channel. The PxPy resonance, which is sit- 
uated just above Ep for Si, is pinned to the Fermi level in 
case of P, therefore the A5 channel becomes highly open. 
For the S impurity the virtual bound state is occupied 
and the transmission of the A5 channel drops. 

To investigate the effect of atomic rearrangements on 
conductance we moved the sp atom from its symmetric 
position. The nanocontact is now simulated by a zig-zag- 
likc chain suspended between Cu leads where the chain is 
a continuation of the fee structure along the [001] direc- 
tion (we show the exact geometry in the insert of Fig. 4). 
The behavior of the conductance vs the atomic mimber 
(Fig. 4) is similar to the previous case: transmission is 
smaller than 1.0 for Na, has a maximum for P and drops 
down for S. However, the absolute values are different 
for the two considered configurations. We found that the 
conductance consists again of three channels. Since the 
symmetry is reduced all channels are different and the 
degeneracy is lost. The 1st and 2nd channels correspond 
to an identity representation (Ag) of the S2 group, while 
the 3rd channel is an antisymmetric one (^„). 



Conductance through 3d atoms attached to Cu 
leads 



As we have seen, when atomic orbitals are present at 
the Fermi level they support additional conducting chan- 
nels. That should be also true for the case of 3d transi- 
tion metal impurities. They have five localized d states 
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purities connected to the Cu leads [see Fig. 1(a)]. Bottom 
panel (b): Symmetry resolved LDOS at impurity sites, at the 
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nearby Ep which potentially might be involved in the 

formation of new channels and the conductance could 
change as compared with « 1.0 Go found for a Cu con- 
striction. Our results for the conductance of the systems 
with 3d transition metal impurities are given in Fig. 5. 
We consider a single atom contact modelled by a short 
atomic chain connecting the Cu leads, with impurity in- 
side the chain (Fig. la). Indeed, we see that conductance 
varies in a broad range along the 3d series reaching val- 
ues of about 2.76*0 for Ti and V, and decreasing to O.SGo 
in the case of a Ni atom. We have excluded many body 
effects (e.g., Kondo effect) from the consideration. As an 
approximation in the calculations the 3d impurities were 
supposed to have zero magnetic moments. That could be 
justified by the fact that at T = the magnetic moment 
of the impurity is effectively screened by the conduction 
electrons. 

Conductance is found to be a sum of three channels 
with Al and A5 symmetry, where the last channel is two 
fold degenerate. To understand the behavior of channels 
transmission (Fig. 5a) we present in Fig. 5b the symmetry 
projected LDOS at the 3d impurities at the Fermi energy, 
as well as in the vicinity of Ep, for selected atoms (Fig. 6). 
Calculations show that the LDOS nearby Ep is formed 
due to 4s and 3d states of transition metal atoms while 4p 
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FIG. 6: Symmetry resolved LDOS at impurity atoms placed 
at the contact site of Cu junction [see Fig. 1(a)]. 



states are above Ep and thus play no role for transport. 

As for 3d states, they form sharp resonances localized in 
the energy window of 1 eV around the Fermi level (Fig. 6, 
low panels). The dxy and dj^2_y2 resonance states, even 
if they are situated at the Fermi level, do not support 
open channels. Since the dxy and d^^-y^ orbitals are 
localized perpendicular to the wire axis effective coupling 
to the neighboring atoms is prevented. Consequently, the 
corresponding resonances look like sharp peaks in the 
LDOS (Fig. 6, down panels). 

Consider first the A5 channel which is projected on 
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dxz, dyz states. The two- fold degenerate resonance state 
of this symmetry is situated above Ep in case of Sc that 
leads to a relatively small transmission of the correspond- 
ing channel. When more d electrons are present at the 
contact atom, the d states move to lower energies. In 
particular, for Ti, V and Cr the dxz,dyz impurity state 
is situated very close to the Fermi level (Fig. 6, down left 
panel) and causes high transmission of the A5 channel. 
For larger atomic number the dxzi dyz resonance moves to 
energies below Ep and transmission of the corresponding 
channel is reduced. 

The highly symmetric Ai channel is open and is mainly 
due to s states for the first half of the 'id series, from Sc 
to Mn. The increase of the s and dz^ LDOS in case of 
Fe, Co and Ni (Fig. 5b) is accompanied by a reduction of 
the transmission. That is, because of the hybridization 
between s states and the d^^ resonance (Fig. 6, down right 
panel) with Ai symmetry, the wave function becomes 
more localized at the impurity atom. Prom another point 
of view, the Ai channel is scattered at the d^^. impurity 
state. For Cu and Zn the d states are shifted to energies 
far below Ep, so that only a single s-channel is left. 



C. Pd constrictions with transition metal 
impurities 

One could expect a variety of different eigenchannels 
for atomic constrictions made of transition metals. An 
example of such a system is a Pd pyramid like contact 
shown in Fig. 1(b), with conductance being 2.8 6*0. The 
LDOS at the contact atom is mainly d-like nearby the 
Fermi energy, as it is seen from Pig. 7. This is also the 
case for the Pd surface. Since d states of different symme- 
try are available at the Fermi energy (Fig. 7), six channels 
give a main contribution to conductance. These are one 
channel of Ai symmetry: ti = 0.72 (spx-like); two dou- 
ble generate channels of A5 symmetry: T2 = T2' = 0.58 
{dxz,dyz-like) and ts = T3' = 0.32 {px,Py-like); and one 
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FIG. 7: Symmetry resolved LDOS at the central atom of the 
Pd junction shown in Fig. 1(b). 



FIG. 8: Conductance through Pd pyramid-like junction con- 
taminated by 3d impurities [see Fig. 1(b)]. 
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FIG. 9: (a) Left column: 
transmission probabilities of 
cigcnchannels through Pd 
junctions contaminated by 3d 
impurities; (b) Right column: 
symmetry resolved LDOS at 
impurity atoms. 



channel of A2' symmetry, = 0.24 (rfj^j^-like). We would 
like to mention, that the channel with the highest sym- 
metry (Ai) has the largest transmission. Here, the irre- 
ducible representations belong to the group. 

Let us now study a situation where the contact atom of 
the otherwise Pd junction is substituted by 3d transition 
metal impurities. The conductance variation along the 
3d series is shown in Fig. 8. We observe that a Ni impu- 
rity does not affect conductance value of 2.8 Go of a Pd 
contact simply because Ni is isovalcnt to Pd. For other 
impurities conductance changes from 1.7 Go in case of Sc 
and Zn, up to 2.9 Go in the case of a Co atom. Eight 
channels give rise to the total transmission. Results of 
the individual transmission probabilities are presented 
in Fig. 9 (left column) together with the LDOS at Ep 
projected on the corresponding orbitals of the impurity 
atoms (right column). Note, that due to a larger open- 
ing angle as compared with the first atomic configuration 
shown in Fig.l, the channels of A2' {dxy) and A2 {d^'^-y'^) 
symmetry contribute significantly to conductance, while 
they were closed in the case of the Cu constrictions with 
2>d impurities discussed in the previous section. 



Without showing a detailed structure of the LDOS 
nearby the Fermi energy, we mention that a high den- 
sity at Ep of the d states of different symmetry is always 
related to resonance-like states (see e.g. Fig. 7 for the Pd 
atom). When some of these resonances approach the 
Fermi level, the transmission of the corresponding chan- 
nel reaches a maximum. Such a situation takes place 
for the A5 channel (d^^z , dy^-like) in case of Ni (Fig. 9, 
middle panels), as well as for the A2 {d^2_y2) and A2' 
(dxy) channels in case of Mn and Co, respectively (Fig. 9, 
lower panels). For pure d channels mentioned above, a 
clear correlation between the partial LDOS at Ep and the 
transmission probabilities is observed. The other channel 
of A5 symmetry does not show a strong variation of the 
transmission vs the atomic number of the impurity. This 
channel is mainly the Px^Py- like. 

Finally, we note that two different channels of Ai sym- 
metry exhibit an interesting variation. One of them is 
spz-like, it is open at the beginning and at the end of 
the 3d series (Fig. 9, upper panels). The dj.2 resonance 
approaching the Fermi level in case of Mn and Fe im- 
purities acts in two ways. Firstly, due to a hybridization 
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FIG. 10: Conductance 
tlirough sp (left plot) and 3d 
(right plot) transition metal 
atoms attached to Cu leads 
in a different way as shown in 
Fig.l. Solid line: pyramid-like 
configuration (Fig. lb), dashed 
line: atomic chain (Fig. la). 
The influence of the electrode 
material on the conductance 
can be seen if we compare 
solid line in plot (b) with 
results of Fig. 8, where a 
geometry of the constriction 
was not changed but Pd was 
used instead of Cu. 



with the s state it reduces the transmission of the already 
existing sp^ channel (that is similar to the case discussed 
in the previous section, Fig. 5). Secondly, the (1^2 reso- 
nance itself opens a new channel of Ai symmetry, which 
consequently is the d-like. Its transmission follows the 
variation of the o?22-LDOS when the atomic number is 
changed. 

D. Does an atom possess certain conductance? 

We have learnt from presented examples that conduc- 
tance properties of individual atoms can be analyzed in 
detail with help of eigenchannel analysis. Therefore the 
question arises: can we attribute definite conductance to 
the particular atom? Comparing transmission through 
3d transition metal impurities attached to Cu and Pd 
leads (Fig. 5a and Fig. 8, respectively), we see that the an- 
swer to this question is in general negative. This state- 
ment is also confirmed by results presented in Fig. 10, 
where we compare transmission through sp and 3c? tran- 
sition metal atoms attached to Cu electrodes in a differ- 
ent way as it is shown in Fig.l. We have mentioned in the 
previous section, that due to a larger opening angle for 
the incident waves, conductance of pyramid-like constric- 
tions increases as compared to nanocontacts simulated by 
atomic chains. The number of non-vanishing conducting 
channels needed to describe the conductance along sp 
and 3c? series becomes also larger. Without going deep 
into details, we point out that number of eigenmodes in- 
creases up to 4 in case of sp impurities (Fig. 10a) and 
up to 6 in case 3c? transition metal atoms (Fig. 10b). In 
particular, in the latter case, because of the increased co- 
ordination number of the contact atom, the d channels of 
(dxy) and A2 {dx2_y2) symmetry become open for Ti 
and Mn, respectively. These channels were closed in the 
case of a nanocontact modelled by a three atom chain 
(see section III.B for discussion). Thus, the presented 
analysis gives no evidence to attribute a precise value of 
conductance to the given atom. A conclusion is that the 



conductance is sensitive to the geometry of atomic con- 
striction as well as to the type of electrodes. However, 
our results show that for the specified single-atomic con- 
tact eigenchannels can always be related to the electronic 
states of the appropriate symmetry at the Fermi energy. 



IV. SUMMARY 

In conclusion, based on the KKR Green's function 
method combined with the Kubo approach we have per- 
formed parameter-free ab initio calculations of the ballis- 
tic conductance through single sp and 3c? transition metal 
atoms attached to Cu and Pd leads. We have investi- 
gated in which way the valency of the atom bridging two 
metallic leads is responsible for the number of conducting 
channels. We have found that both the chemical valency 
of the atom and the opening angle of the atomic constric- 
tion determine the number of eigenmodes. In accordance 
with the known tight-binding model by Cuevas et al.^ we 
have confirmed by our study that the symmetry of the 
open conducting channels is related to the symmetry of 
the electronic states at the contact site. That allows to 
relate the open channels to the orbitals of the contact 
atom available at the Fermi energy. Furthermore, we 
have shown that impurity resonances approaching the 
Fermi level can open new channels of the appropriate 
symmetry. On the other hand, in some cases, the trans- 
mission of the existing channels can even be blocked by 
scattering at virtual bound states which have in this case 
the symmetry compatible with the symmetry of the open 
channels. 
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